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O . We investigate the scattering off three nonoverlapping disks equidistantly spaced 

along a line in the two-dimensional plane with the radii of the outer disks equal 
^ . and the radius of the inner disk varied. This system is a two-dimensional scattering 

analog to the double-well-potential (bound state) problem in one dimension. In 
both systems the symmetry splittings between symmetric and antisymmetric states 
or resonances, respectively, have to be traced back to tunneling effects, as semiclas- 
sically the geometrical periodic orbits have no contact with the vertical symmetry 
axis. We construct the leading semiclassical "creeping" orbits that are responsible 
for the symmetry splitting of the resonances in this system. The collinear three- 
disk-system is not only one of the simplest but also one of the most effective systems 
for detecting creeping phenomena. While in symmetrically placed ra-disk systems 
creeping corrections affect the subleading resonances, they here alone determine 
the symmetry splitting of the 3-disk resonances in the semiclassical calculation. It 
should therefore be considered as a paradigm for the study of creeping effects. 
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1 Introduction 



Quantum mechanics tells us that the eigenf unctions of the one-dimensional Schrddinger 
equation of lowest eigenvalue have no nodes, the next-lowest eigenfunctions have one 
node, etc. So, it is well known that the ground state of a particle in the (symmetric) 

V(x) 




Figure 1: Symmetric double- well potential V(x) plotted (solid line) together with the corre- 
sponding single- well subpotentials (dashed lines). In addition, the classical turning points — X2, 
—x\, x\ and X2 at the total energy E are shown and the geometrical periodic orbits (long-dashed 
line) as well as the tunneling orbit (short-dashed line) are indicated. 

double- well potential (see Fig.[l|) is not centered at the bottom of one well, as one might 
naively infer from purely classical considerations, but instead is given by the spatially even 
combination of (approximately) harmonic-oscillator states centered at the bottom of the 
two wells, whereas the first excited state is given by the spatially odd combination. The 
degeneracy of the two ground-state energy eigenvalues belonging to each of the single- well 
potentials individually is broken in the double-well case. However, this splitting cannot 
result from any perturbative h corrections that alter the even and the odd combinations 
in the same way, but only from the barrier penetration. In fact, the difference between 
the odd and even energy combination is proportional to the barrier-penetration factor 

where —X\, x\ are the (inner) classical turning points, V(x) the potential and E the 
ground-state energy, and m is the mass of the particle. 

The splitting is thus inherently a nonperturbative effect linked semiclassically to tun- 
neling orbits or instantons; see, e.g., Refs. Jj], 0, [3|. This means that periodic orbit theory 
with only geometrical orbits cannot describe the splitting of the single-well states in even 
and odd double-well states, as long as the energy is below the potential barrier. Periodic 
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orbits with tunneling sections have to be added to the semiclassical theory. The reason 
is that as long as the energy is below the potential barrier the geometrical orbits can 
never reach the symmetry axis of the double-well potential at x = and therefore are not 
sensitive to the boundary condition (Neumann or Dirichlet) chosen there. 

Here we want to construct the simplest two-dimensional and scattering analog of the 
double-well-potential problem. It consists of three nonoverlapping disks centered and 
equally spaced on a straight line in the two-dimensional plane, where the outer disks have 
the same radius, whereas the radius of the inner disk is free to vary [see Fig. 0(a)] . 




Figure 2: (a) Collinear three-disk system in the full domain. The outer disks have a radius 
of size a, whereas b is the radius of the inner disk. The center-to-center separation R and the 
length 2L of the geometrical periodic orbit (thin solid line) are related as L = R — a — b. (b) 
The corresponding two-disk (sub)system in the full domain, (c) The half-symmetry reduced 
collinear three-disk system in the half domain. In addition, the vertical symmetry axis is shown 
and the leading creeping periodic orbits are indicated by dotted and dashed lines. 

The analog of the single-well potential is the two-disk subset of the upper problem 
consisting of one of the original outer disks and the middle disk, where the radii as well 
as the distances of the disks are kept unaltered [see Fig. 0(b)]. As in the double- well 
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potential, the geometrical periodic orbits never hit the vertical symmetry axis, which 
goes through the center of the middle disk perpendicular to the horizontal symmetry 
axis (the latter is in turn given by the line through the centers of the three disks); see 
Fig.^|(c). The splitting can only come from nongeometrical periodic orbits. In fact, only 
creeping periodic orbits are left over for this purpose (see Refs. 0, ||, [7], §], based on 
Refs. |9|, [10]] ) . These simple phenomena should not be mixed up with the chaos- assisted 
tunneling of Refs. [11, [12[] or the geometrical interpretation of multidimensional tunneling 
for bound state systems of Refs. O, [TJ . 



In Sec. 2 we will describe the collinear three-disk system in more detail and in Sec. 3 
we outline the calculation (details can be found in the Appendix) and give results. Con- 
clusions are given in Sec. 4. 



2 Collinear three-disk system 

The standard three-disk system with its three equally sized and nonoverlapping disks cen- 
tered at the corners of an equilateral triangle in the two-dimensional plane is discussed in 
the literature [|TJ, [16|, [L7|, [18], [Lj], |20], |21] as one of the simplest (if not the simplest) classi- 
cally completely chaotic scattering system. Here we consider a different three-disk system 
with the centers of the three disks arranged in equal intervals on a straight line, where 
the outer disks have equal radii a, whereas the radius b of the inner disk is free to vary; 
see Fig. 0(a). As the original three-disk system is the paradigm for a chaotic scattering 
system, the new collinear three-disk system can be taken as the paradigm for scattering 
system with maximal creeping corrections. Note that the collinear three-disk system is 
invariant under reflections with respect to the horizontal symmetry axis (defined by the 
centers of the three disks) and - as long as the outer disks are of the same radius and the 
spacings between the disks are equal - with respect to the vertical symmetry axis (which 
goes through the center of the middle disk, perpendicular to the horizontal symmetry 
axis). Thus the states and resonances of the collinear three-disk system can be classified 
according to the four irreducible representations of the the Civ group, whereas each of the 
two-disk subsystems in general only has a ah reflection symmetry with respect to its hori- 
zontal symmetry axis (as the radii of the two disks can be different). As usual, it is easier 
to study the quantum mechanics and semiclassics of these systems desymmetrized in the 
fundamental domain, instead of the original systems in the full region. The fundamental 
region of the collinear three-disk system is only one quarter of the full region, bounded 
by the horizontal and vertical symmetry axis - see Fig. 0(b). The totally symmetric A\ 
(^-representation is characterized by even or Neumann boundary conditions for the wave 
functions on the horizontal and vertical symmetry axis, whereas the totally antisymmet- 
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horizontal symmetry axis: Al: Neumann boundary conditions 

A2: Dirichlet boundary conditions 




vertical symmetry axis: 

Al & Bl: Neumann boundary conditions 

A2 & B2: Dirichlet boundary conditions 



horizontal symmetry axis: Al & B2: Neumann boundary conditions 

A2 & B 1 : Dirichlet boundary conditions 

Figure 3: (a) Two-disk (sub)system in its fundamental domain, (b) Collinear three-disk system 
in its fundamental domain. Indicated are the geometrical periodic orbits (solid lines) and the 
leading creeping periodic orbits (dashed and dotted lines). In addition, the irreducible C 2 
and C 2v representations of the two- and three-disk system, respectively, and the corresponding 
boundary conditions on the horizontal and vertical symmetry axes are shown. 

ric A 2 representation has just the opposite transformation behavior, i.e., the boundary 
conditions are both of the Dirichlet type. In addition, there exist two more irreducible 
representations of mixed symmetry, the Bi and B 2 representation; the former transforms 
symmetrically with respect to the vertical symmetry axis and antisymmetrically with the 
respect to the horizontal one and the B 2 representation has the opposite transformation 
properties. The two-disk subsystems have only two irreducible C\h (— C\ v = C s ) repre- 
sentations: the totally symmetric A' one, which is even with respect to reflections off the 
horizontal symmetry axis, whereas the states of the totally antisymmetric A" represen- 
tation have an "odd" transformation law. In the following we will use the nonstandard 
notation A\ (A 2 ) instead of the standard on A' (A") for the C\h representations; see, e.g., 



Ref . p2f . The fundamental region of the two-disk subset is therefore the half plane with 
respect to the horizontal symmetry axis; see Fig.|3](a). 

Note that the collinear three-disk system possesses only two geometrical periodic orbits 
in the full domain. They are symmetric with respect to the vertical symmetry axis and 
therefore map into only one geometrical periodic orbit in the fundamental domain [see 
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Figs. 0(a) and 0(b)]. It runs along the horizontal symmetry axis and should therefore 
be treated as a boundary orbit (see Refs. p3[ for a discussion of boundary orbits). 
The same is true for the solely existing boundary orbit of the two-disk subsystem, which 
also runs along the horizontal symmetry axis. This means that the contributions of the 
geometrical orbit to the B\ and A 2 states and resonances of the collinear three-disk system 
and to the A 2 resonances and states of the two-disk subsystems are strongly suppressed 
[see Eq. (|A7|) of the Appendix] . 

As the symmetries with respect to the horizontal symmetry axis are common for the 
collinear three-disk system and its two-disk subsystems, they are not essential for our 
discussion. Instead of working with the full three-disk-system, one could as well only 
consider it in the half domain defined by its horizontal symmetry axis [not by its vertical 
one as in Fig. 0(c)], and compare it to the two-disk subsystem in the fundamental domain 
[Fig.|3|(a)]. By choosing either Neumann or Dirichlet boundary conditions on the hori- 
zontal symmetry axis of these system one can then enhance or suppress the geometrical 
relative to the creeping contributions. So in an experimental billiard setup, one could 
work with half disks (instead of full disks) that are aligned at the reflecting boundary of 
an otherwise "open" billiard geometry; see, e.g., Ref. [25|. To our knowledge this reference 
contains the only known experiment on the scalar open n-disk scattering problem. It is 
done with microwaves in a cavity with absorbing outside walls. As long as the extension 
in the third direction (perpendicular to the quasi-two-dimensional n-disk setup that is 
experimentally realized with reflecting metal cylinders) is so small that only the lowest 
mode is excited in that direction, the stationary electromagnetic problem is mathemat- 
ically equivalent to the stationary two-dimensional quantum (scattering) problem of a 
point-particle. In that respect this experiment as well as any possible extension to the 
collinear three-disk setup serves just as an "analog computer" to our digitally computed 
data. 

The analog of the original one-dimensional double- well potential is the partially desym- 
metrized collinear three-disk system in its half domain relative to its horizontal symmetry 
axis and the analog of the two single-well potentials are the two two-disk subsystems 
in the fundamental region. The single- well eigenfunctions correspond to the A\ (or A 2 ) 
states of the two-disk subsystems, and the even and odd combinations of the single-well 
functions to the A\ and B 2 (or B\ and A 2 ) states of the collinear three-disk system in 
its fundamental domain [Fig.^(b)]. The correspondence of the geometrical orbits is ob- 
vious. In the double well as in the collinear three-disk system the geometrical orbits do 
not reach the "x=0 axis" or the corresponding vertical symmetry axis, respectively. They 
cannot produce the symmetry splitting. The role of the tunneling orbit in the double-well 
potential is taken over by those creeping orbits of the three-disk system that reach the 
vertical symmetry axis [see, e.g., Figs.^(c) and^(b)] as only those orbits can produce the 
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symmetry splitting of the A\ {A2) two-disk resonances into the A\ and B2 {B\ and A2) 
three-disk resonances, corresponding to the splitting of the single-well eigenenergies into 
the even and odd double-well eigenenergies. However, whereas in the double-well problem 
the tunneling corrections scale exponentially with h~ , the creeping corrections along a 
circular path scale exponentially with h~z |], |10[| . Common to both is the nonperturba- 
tive structure. In detail, the creeping tunneling exponent of mode number £ = 1, 2, 3, . . . 
is given as [JTOj, || 



Here p is the momentum, p(s) is the local radius of curvature of the creeping path (which 
parametrically depends on the length s of the creeping arc), and qe (with t = 1, 2, 3, . . .) 
is the £ th zero of the Airy integral A(q) = J °° dt cos(tg — t 3 ) which can be approximated 

as 

r 



qe ~ -63 



3tt 



(3) 



If x a and Xh are the start and end points of the creeping section, the semiclassical creeping 
Green's function reads 

G?(x a ,x b ,p) = e -Io ds '^ s '^e^ Xa ' Xb ^ , (4) 

where s is the length of the creeping arc from the point x a to the point Xb and S(x a , Xb,p) 
is the action along it. For the special case that the creeping arc is of circular shape of 
radius a with an arc angle 5<j), the formula (||) then simplifies to 



Gf(6<p,p) = exp (- gi e-^ 6 (J£) 3 6<j>) exp QW) , (5) 



such that the nonperturbative creeping contribution can be read off easily. Note that the 
exponent of the (leading) creeping "tunneling" suppression factor exp[—gi cos(7r/6)(/ca/6) 36(f)} 
(where k is the wave number) scales linearly with the creeping angle and only as third 
root of the creeping radius. Thus the creeping suppression is governed by the creeping 
angle and the size of the disk is only of secondary importance. We will observe this fact 
in the comparison of the exact to the semiclassical data (see also the discussion at the 
end of the Appendix). Note further that the creeping is suppressed with increasing wave 
number k, such that the splitting eventually vanishes in the semiclassical limit Re k — » 00. 



3 Calculation and results 

In this section the quantum-mechanical and semiclassical calculations are briefly explained 
(details can be found in the Appendix) and finally results for the splitting of the two-disk 
resonances into the resonances of the collinear three-disk system are shown. 
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As the collinear three-disk system and the two-disk subsystems involve disks of dif- 
ferent sizes, the quantum-mechanical calculation of Ref. NTH] has to be generalized. This 



was done in Ref. [26], where in fact the characteristic determinant detNl(k) of any n-disk 
system involving n < oo nonoverlapping disks of any sizes has been constructed (see the 
Appendix). The zeros in the lower complex wave-number (k) plane of this Korringa- 
Kohn-Rostoker-type f27 , |28|| determinant define the genuine multidisk scattering reso- 
nances [18, 26 1 . The fact that the collinear three-disk system and its subsystems have 
disks of two different sizes implies that the corresponding characteristic matrices M(A;), 
even in its desymmetrized form, cannot just be expanded in the angular momentum 
eigenbasis {| I)} of one disk alone, but have to be expanded in the so-called superbasis 
{|i;0}> which acts on the disk surface sections j in the fundamental domain [f26] 
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also the Appendix). Thus the cumulant expansion of the characteristic determinants (see 
Refs. @, [|, ^6|) is effectively not organized in terms of the first cumulant TrA(A;) [where 
A(k) = M(A;) — 1 is the trace-class kernel [PUfl l, but in terms of the second and higher 
cumulants, e.g., |{(TrA(A;)) 2 — TrA 2 (A;)}, etc. The semiclassical contribution of the first 
cumulant in the desymmetrized collinear three-disk system in fact corresponds to a ghost 
orbit (see, e.g., Refs. [^9], pB[ ) between the outer two disks, which is obstructed by the 
presence of the middle disk. All the "nonghost" geometrical and creeping periodic orbits 
of the collinear three-disk system and its two-disk subsystems have therefore an even topo- 
logical length and especially the shortest orbits in the fundamental domain are already of 
topological length 2. From the semiclassical point of view of Refs. [|], [5| this is obvious as 
the orbits have an even number of contacts (in the sense of specular reflections for geomet- 
rical legs and tangential creeping for creeping legs) with the disks. The number of these 
contacts corresponds just to the power m in the trace TrA m (/c), as the m th trace involves 
m sums. In the semiclassical calculation, each of these m sums are replaced either by an 
integration and a following saddle point approximation (which corresponds to a specular 
reflection) or by a sum over the creeping poles (which corresponds to a creeping contact). 

As with increasing topological order the number of creeping orbits is increasing dra- 
matically and as, on the other hand, the leading creeping orbits give the dominant contri- 
bution to the splitting of the resonances, we show only semiclassical results for the lowest 



topological order, namely order 2. This corresponds to the curvature expansion [30, pT[ 



of the Gutzwiller-Voros zeta function [^, up to this topological order [see Eqs. ( |A14| ) 
and (|A15|) of the Appendix]. For a qualitative comparison to the exact data this is more 
than sufficient. 

In Figs. f| and [5] the exact collinear three-disk A\ and B 2 resonances (plotted both 
as diamonds) are compared with the exact A\ resonances of the two-disk subsystem 
(plotted as crosses). The splitting in the imaginary part is clear from the figures. The 
splitting in the real part is of similar size, but compressed in the figures because of the 
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different scales used for the real and imaginary &;-axes. In addition, the predictions of 
the geometrical orbit [see ( P^7| ) and ( |A15| )] are plotted (as boxes), which describes very 
well the leading resonances of the two-disk subsystem, with more and more accuracy for 
larger and larger Re k. The deviations between the exact quantum-mechanical data and 
the semiclassical predictions which are visible in Figs.|l-|8] for small wave numbers result, 



in the case Rek > 1/a, from the neglect of % corrections 0, [35], ^6|, subleading creeping 
orbits (see the discussion at the end of the Appendix) and subleading terms in the Airy 
expansion of the creeping propagators |37|, |38| . The deviations for Re A; <C 1/a can be 
traced back to a breakdown of the standard semiclassical expansion in terms of geometrical 
and creeping orbits altogether. In that wave-number regime the quantum-mechanical 
wave cannot resolve the scattering obstacles, so that the semiclassical expansion has to be 



expressed in terms of diffractional orbits scattered from point-like centers; see Ref. [39 
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Figure 4: Exact quantum-mechanical (A\ and B2) resonances of the collinear three-disk system 
(with b = 1.5a and L = 4a) shown as diamonds in the complex wave-number (k) plane and the 
A\ resonances of the corresponding two-disk subsystem are denoted by crosses. The predictions 
of the geometrical orbit [Eqs.([A7|) and ( A15| )] are presented by boxes. The real and imaginary 
parts of k are measured in units of 1/a. Note the different scales for the real and imaginary k 
axes. 
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In Figs. ^ and [5] the A\ and B 2 resonances of the collinear three-disk system are compared 
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Figure 5: Same as in Fig.^, however, with b = 3a. 



with the semiclassical calculation ( |A14| ) based on the geometrical orbit ( |A7| ) and leading 
creeping corrections ( [All]) and (|A12 ). The creeping terms reproduce, at least qualita- 
tively, the trend of the data, i.e., the tunneling splitting between the A\ and B 2 resonances 
of leading order. The subleading A\ and B 2 resonances are of course not reproduced, as 
they are determined by higher-order contributions (geometrical plus creeping ones) in the 
curvature expansion [JTJ of the semiclassical Gutzwiller-Voros zeta function 
Furthermore, the stronger deviations for the low lying resonances (i.e., the systematical 
underestimation of the magnitude of the imaginary part) can be explained by the ne- 
glect of higher correction terms (beyond those discussed in Ref. in the Airy expansion 



of the creeping propagators; see, e.g., Ref. [3/] for Airy correction terms in the one-disk 



scattering system and [38 for results in the two-disk system. 



Finally, in Fig. |] the exact quantum-mechanical A2 and B\ resonances of the collinear 
three-disk system are compared with the corresponding semiclassical predictions calcu- 
lated from Eq. ( |A 1 4j ) . Furthermore, the exact A 2 resonances of the two-disk subsystems 
are shown together with the semiclassical prediction from Eq.( |A15|) , which involves only 
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Figure 6: Exact quantum-mechanical (A± and B 2 ) resonances of the collinear three-disk system 
(with b = 1.5a and L = 4a) shown as diamonds and boxes, respectively, in the complex wave- 
number (k) plane. The corresponding semiclassical results of Eq.(A44), which include the two 
leading creeping orbits, are shown as upright and diagonal crosses, respectively. The real and 
imaginary parts of k are measured in units of 1/a. Note the different scales for the real and 
imaginary k axes. 

the geometrical orbit. By comparing these resonances with Figs.^ and ^ (which describe 
the A\ and B 2 three-disk and A\ two-disk resonances) one can learn that the resonances 
of Fig.H are suppressed. This result is expected, as the geometrical orbit now runs on 
a Dirichlet line that bounds the fundamental domain (of the collinear three-disk as well 
as the two-disk system), such that the wave function goes to zero there, whereas before 
the geometrical orbit was affected by the Neumann condition. Note that the two-disk 
resonances are again quite well approximated (with increasing Re k) by the predictions of 
the geometrical orbit alone. (The agreement is worse than in Figs. ^| and [7] as the reso- 
nances here are stronger suppressed to begin with, such that neglected subleading effects 
are of relative higher importance. In fact, most of the deviations should be traced back to 
the neglected creeping orbits of the two-disk system; see, e.g., Ref. [JJ for the discussion 
of the similar B\ resonances of the symmetric two-disk system.) However, the inclusion 
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Figure 7: Same as in Fig.|| however, with b = 3a. 



of the creeping orbits in the semiclassical calculation is necessary in order to predict the 
qualitative trend of the exact three-disk data. Furthermore, one can learn from these data 
that the creeping orbits are now not only leading in the splitting, but in fact they even 
dominate the geometrical orbit in the determination of the absolute position of the reso- 
nances. This is, of course, again a consequence of the Dirichlet choice for the boundary 
condition on the horizontal symmetry axis, which suppresses the geometrical contribu- 
tions relative to the creeping ones. In the case where one wants to maximize the creeping 
effects, the Dirichlet choice for the boundary condition on the horizontal symmetry axis 
is advantageous. However, one should note that this choice suppresses the resonances 
altogether, such that the Neumann choice might be still better from the phenomenolog- 
ical point of view, as the resonances can be more easily identified. By positioning two 
receiving antennas symmetrically to the horizontal symmetry axis the even and odd states 
with respect to that axis can be extracted experimentally. The "antennas" refer of course 
to the electromagnetic two-dimensional analog case (see, e.g., Ref. p5| ) of the quantum 
scattering problem in two dimensions. By adding another pair of antennas symmetrically 
to the vertical symmetry axis, the remaining symmetries can be determined. 
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Figure 8: Exact quantum-mechanical (A2 and B\) resonances of the collinear three-disk system 
(with b = 1.5a and L = 4a) shown as diamonds and boxes, respectively, in the complex wave- 
number (k) plane. The corresponding semiclassical results of Eq.(A14), which include the two 
leading creeping orbits, are shown as upright and diagonal crosses, respectively. Furthermore, 
the exact A2 resonances of the two-disk subsystem are plotted as triangles and the corresponding 
semiclassical predictions based solely on the geometrical orbit [see Eq.( A15| )j as stars. The real 
and imaginary parts of k are measured in units of 1/a. Note the different scales for the real and 
imaginary k axes. 



4 Conclusions 



We have shown that the collinear three-disk system is a two-dimensional scattering analog 
(probably one of the simplest) of the one-dimensional two-potential bound-state problem. 
As in the latter case, nonperturbative contributions are needed in order to describe the 
splitting of the resonances (states) of the subsystems (the two-disk or one-potential prob- 
lems, respectively) into the corresponding quantities of the full system. Whereas in the 
two-potential problem these nonperturbative contributions are given semiclassically by 
tunneling paths, in the collinear three-disk scattering system their role is taken over by 
creeping orbits, which alone determine the splitting phenomena. In the description of the 
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leading resonances of the collinear three-disk system they are at least as important as the 
geometrical orbit. In the case where the Dirichlet boundary condition has been specified 
on the horizontal symmetry axis the creeping orbits even dominate the geometrical or- 
bit. Thus the creeping orbits of the collinear three-disk system are of "leading nature". 
This is the main qualitative difference to the standard two- and three-disk systems with 
symmetrically arranged disks, as there the creeping terms are inessential for the descrip- 
tion of the leading band of resonances and do only play an important role in the case of 
subleading bands. In summary, the collinear three-disk system is probably the simplest 
(higher than one-dimensional) scattering analog of the two-potential bound-state prob- 
lem and the simplest disk-scattering problem with leading creeping terms. In the same 
way as the standard two-disk problem is the paradigm for hyperbolic scattering systems 
and the standard equitriangular three-disk system is the paradigm for chaotic scattering 
systems, the collinear three-disk system should be considered as a paradigm for the study 
of creeping effects. 

By replacing the disks in the collinear three-disk system by three-dimensional balls, 
one can easily construct a three-dimensional scattering analog system to the double-well 
potential problem. As in the two-dimensional collinear three-disk case, the splitting of 
the resonances (with respect to the vertical symmetry plane) again has to be traced back 
to the creeping orbits. 

As mentioned in Sec. 3, the geometrical orbit spanned by the outer two disks is shad- 



owed by the presence of the middle disk and is therefore only a ghost orbit [^, |28[ . This 
ghost orbit can be tested in a modified three-disk array, where the middle disk is posi- 
tioned off the horizontal symmetry axis along the vertical one. If the middle disk does not 
overlap the old center of the collinear three-disk system, the ghost orbit in fact splits into 
the sum of a direct geometrical orbit that is spanned by the outer disks and a geometrical 
orbit that, in addition, encounters a specular reflection from the middle disk. By sliding 
the middle disk back into its old position, one can check how and when exactly the two 
geometrical orbits become a ghost orbit (see, e.g., Ref . [|4~0fl for the introduction of glanc- 
ing orbits [fll]] into periodic orbit theory). This phenomenon is strongest if the middle 
disk possesses Neumann boundary conditions, as then the two geometrical orbits inter- 
fere constructively, whereas they interfere destructively in the case of Dirichlet boundary 
conditions. In this respect a very interesting system would be a nonoverlapping four-disk 
system with the centers of three equally sized disks at the corners of an equilateral trian- 
gle, whereas the fourth disk (again with Neumann boundary conditions) is positioned at 
the center of this triangle. If the size of this fourth disk is varied, geometrical orbits can 
be turned into ghost orbits or vice versa. Furthermore, if the fourth disk is large enough, 
again only creeping orbits are left over for producing the splitting of the resonances of the 
three two-disk subsystems (consisting of one outer disk and the middle disk and classified 



14 



according to the point-symmetry group C\h) into the resonances of this new four-disk 
system, which are classified according to the point-symmetry group C 3v . Whereas the 
collinear three-disk system is classically hyperbolic, but nonchaotic (as its two-disk sub- 
systems), the new four-disk system is classically completely hyperbolic and chaotic if the 
middle disk is small enough. By varying the ratio of the size of this fourth disk relative to 
the common size of the other three, one has the possibility to study the transition from 
an nonchaotic to a chaotic system classically, semiclassically and quantum mechanically. 
In the case where the boundary conditions on the fourth disk are also of Dirichlet nature, 
the formulas of Ref. suffice to attack the quantum mechanics of the problem. Again 
three-dimensional generalizations are easy to generate, e.g., the four-ball analog system 
or, even more interesting, a (nonoverlapping) five-ball system in tetrahedral form, where 
four balls are positioned at the corners and the fifth in the middle. 

All these suggested systems can also serve as a higher- dimensional scattering analog 
of the double- well potential. Still, the collinear three-disk system is by far the simplest. 
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A Appendix 



The characteristic n-disk matrix M of Ref. [26| reads, in the full n-disk domain, 



M# = 8»'6 W + (l - S»') ^ MW^'-^'H-l) 1 ', (Al) 



where — oo < 1,1' < +00 are the angular momentum quantum numbers in two dimensions, 
j,j' label the disks, the J;(fcr)'s are ordinary Bessel functions, and the H^^krYs the 
corresponding Hankel functions of first kind. The quantity aj is the radius of disk j, Rjj> 
is the distance between the centers of disks j and j', and atfj is the angle of a ray from the 
center of disk j to the center of disk j' as measured in the local (body-fixed) coordinate 
system of disk j. These quantities take the following values for the collinear three-disk 
systems: j = 1, 2 labels the outer disks on the right and left side, respectively, whereas 
j = is reserved for the inner disk. Furthermore, we have 

a 2 = a\ = a , 
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-Roi — -Rio — R20 — Rq2 — R 1 

R12 = R21 — 2-R , 

«21 = «12 = , 

«20 = «02 = 7T , 

o;oi — ft but «io = 0. 



(A2) 



The characteristic matrix M of the two-disk (sub)system is of course a special case of the 
above, namely either j — 0, 2 or j — 0, 1. 

The corresponding desymmetrized M matrix of the collinear three-disk system in the 
fundamental domain reads 

2 



M°?[7J] 

MM 
MM] 



l + (_l)/+^ /! + (_!)/+!+«' 



l + ^o 
V2 



l+(-l) I+1+l b Ji(kb) 



\I+J+l' o-(l) 



H\+ v (kR) 



(-l)-'v^ 'a J,(fca) ^(.^^ 



+ 



;i+^o)(i+^ ) 

l+(-l)^o 
1+Sio 



Ji(ka) 



\l+S l0 )(l+S ll0 )H^(ka) 



bH\}\kb) 



Hi%{2kR)H-l) I+J+l ' H^kR)] , (A3) 



with </,/'< 00 and where the indices I and J specify the irreducible Ci v representations 
with the following identification: I — 1,2 for the Aj, Bj representations, respectively; 
i.e., (— 1) J+1 is the phase under reflection off the vertical symmetry axis and (— is 
the phase under point reflection with respect to the center of the 3-disk system. 

The characteristic matrix of the desymmetrized two-disk subsystem in the fundamental 
domain reads 



MM 
MM 

mm 
mm 



5u> 



1 + 



-l) K+ %o 



1 



l + <S,o)(l + <W a H^\ka) 
{-l) l+v a J t (ka) 



b Mkb) [Hn(kR) + (-lf^'H^(kR) 
H\%{kR) + (-l) K +^'Hll],(kR) 



1 + (-l)* +1 ^o 



1 + 5, 



10 



(A4) 
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with <l,l'< oo and where the index K — 1,2 refers to the irreducible Ak representation 
of the point-symmetry group Cih- More specifically, (— 1) K+1 = (— 1) I+J is the phase 
of the two-disk and three-disk states under reflection off the horizontal symmetry axis. 
Finally, the resonances k res are determined as the zeros of the desymmetrized characteristic 
determinant [evaluated in the superspace, {|j, /)} with j=0,l and < / < oo, acting on 
the disk surfaces (see Figs. U) of the fundamental domain] in the (lower) complex wave- 
number (k) plane 

{detM[/jp)| _ =0 (A5) 

in the collinear three-disk case and 

\detM[K](k)\ _ =0 (A6) 

in the case of the two-disk subsystems. 

The fundamental geometrical orbit that is of topological length 2 and that is common 
to the collinear three-disk system as well as to its two-disk subsystem is given by 

i2Lk 

t G (k; K) = -== -. , (A7) 

where the index K — 1,2 has been defined above. Thus the K = 2 representations, 
namely, the A 2 two-disk representation and the A 2 and B\ three-disk representations, are 
suppressed relatively to the K = 1 representations, namely, the Ai two-disk one and the 
Ai and B 2 three-disk representations. This is intuitively clear as the geometrical orbit is 
a boundary orbit on the horizontal symmetry axis and therefore sensitive to the choice of 
the boundary condition for the wave function, i.e., either Dirichlet or Neumann boundary 



conditions; see Refs. |23|, |24}| . The length 2L of the geometrical orbit is given by 

2L = 2(R-a-b) , (A8) 

where a and b are the radii of the outer and inner disk, respectively, and R is the center- 
to-center separation between these disks. The quantity A is the leading eigenvalue 

/ / ab\ 2RL . K . 

A = 1+ ( 1 + V 1 + 7irJ^r < A9 > 

of the stability matrix, which in turn is given by the product of four submatrices: a 
translational one from the inner to the outer disk, then a reflectional matrix linked to the 
outer disk, then again a translational one from the outer to the inner disk, and finally 
again a reflectional matrix acting at the inner disk 

^ I I ! ! I ! s / f I I S I II i f I '/ft I 

(AlO) 





17 



(see Ref. [^T] for further details). Note that the reflection angle 6 is zero for the two reflec- 
tions, such that the off-diagonal elements 2/[pjCos(#)] of the reflection matrices simplify 
to 2/ pi where pi is the local radius of curvature. As there are two reflections and as the 
geometrical orbit does not touch the vertical symmetry axis, the Maslov phase as well as 
the group-theoretical phase are zero (modulo 2tt) in (|A7|) . 

Under the so-called Keller construction of Ref. [0] the two leading creeping orbits of the 
collinear three-disk system in the fundamental domain [see Fig. 0(b) , dotted and dashed 
lines] have the structure 

/ h p i7T / 12 r ,pi{k2Ll+54>i[kb+q t e™/*{kb/&)l)} 

t cl (k;J,K) = -(-l) J +\ -^- r y±£^ r> (All) 



tc2(k;J,K) = -(-^) J+K \l^ L ^ < —^Y — — — — - (A12) 

v ; y ! V 2i?f (kb) l e ti i - ^-iy+K el n[kb +qe e^ { kbmh v ; 

The creeping parameters qg and Cg are defined in Ref. 0, see also Ref. qi is the £ th 
zero of the Airy integral, A(q) = / °° dt cos(qt — t 3 ), and 

w — — — 



m/12 Q e i{k2Ll+8<t> 2 [k b +q e e in ^(kb/6)^]} 



The indices J and K have been defined above. Furthermore, there enter the lengths Lf 
of the geometrical legs, the creeping angles 5<pi, and the effective radii R^ s of the creeping 
orbits (in the notation of Ref. ||). The former two quantities as well as the specular 
reflection angles 9i can directly be read off from the geometry of Fig. 0(b). The effective 
radii are constructed by utilizing the formula (see Ref. ||) 

m 

R eS = £oH(l + £nK n ) , (A13) 

n=l 

where £ n is length of the geometrical leg between the n th and (n + l) th points of reflection 
and n n is the curvature right after the n th reflection. Here m = 1, £o(i) = £i(i) = Lf and 
Ki(i) = (1/Lf) + 2/ (a cos #,). In summary, we have the following expressions: 

Lf = 

f 7i I b \ I 

< 71 



0i 




< 71 
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cff 




L g 

•> i - ii 

\ acos{U 2 ) 



Rf = 2L\ 1 + 



Note the minus signs on the right-hand side of ( |A11| ) and ( |A12| ), as both orbits encounter 



only one specular reflection off a disk in the fundamental domain. The phases (— 1) J+1 
and (— 1) J+K are responsible for the splitting of the Ai and B 2 (Bi and A 2 ) three-disk 
states relative to the A\ (A 2 ) two-disk states, respectively, and result from the number of 
contacts of the creeping orbits with the vertical and horizontal symmetry axis, respectively. 

Equations ( |A1 1| ) and (|A12|) are the only creeping orbits in the fundamental domain of 
the three-disk (and also two-disk) system that have a potential creeping angle 5<fri <C tc (in 
case ici?). The other creeping orbits of the two- and three-disk system have creeping 
angles of at least 5(f>i > n and are therefore strongly suppressed relative to the above 
two. In fact, in Ref. M it was shown for the standard two-disk system that periodic 
orbits with creeping sections (which all have creeping angles 5<pi > it) hardly affect the 
leading resonances on which we concentrate here. They do give, however, appreciable 
corrections to the nonleading resonances. So as long as we limit our discussion to the 
leading resonances, we only have to take into account the creeping orbits ( |A11|) and ( |A12|) 
in addition to the geometrical orbit (|A7|) , of course. The inclusion of the left-out creeping 
orbits (of topological length 2) in the semiclassical calculation shifts only the first two 
leading resonances of Sec. 3, whereas all the other leading resonances are left unchanged 
up to figure accuracy. For the same reasons, we can also neglect the I > 1 creeping modes 
and the creeping terms in the denominator of ( [All] ) and ( |A12| ), respectively. Thus the 
leading semiclassical resonances are determined from the following relations, which are 
the first curvature approximations (they actually corresponds here to topological length 2, 
see Sec. 3) to the Gutzwiller-Voros zeta function with and without diffraction corrections, 
respectively (see, e.g., Refs. [|, p6| ) 

= {l-t G (k;K)- t cl (k; J, K) - t C2 (k- J, K)} k=k _ (A14) 

for the At (J=l, K=l), A 2 (J=2, K=2), B x (J=1,K=2) and B 2 (J=2, K=l) resonances 
of the collinear three-disk system and just 

= {l-t G (k;K)} k=kree (A15) 

for the Ax resonances (K—1,2) of the two-disk subsystem. 
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